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Abstract 

We study three-dimensional deformations of thin inextensible elastic rods witli non- 
vanishing spontaneous curvature and torsion. In addition to the usual description in 
terms of curvature and torsion which considers only the configuration of the centerline 
of the rod, we allow deformations that involve the rotation of the rod's cross-section 
around its centerline. We derive new expressions for the mechanical energy and for 
the force and moment balance conditions for the equilibrium of a rod under the action 
of arbitrary external loads. Several illustrative examples are studied and the con- 
nection between our results and recent experiments on the stretching of supercoiled 
DNA molecules is discussed. 

87.15.La, 46.70.Hg 
1. Introduction 

Recent experimental advances in the art of manipulation of single DNA molecules and of 
rigid protein assemblies such as actin filaments, etc., have led to an outbreak of theoretical 
activity connected with the elasticity of thin rods^~^^. One of the most intriguing theoretical 
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questions related to the deformation of DNA concerns the couphng between bending and 
twist in the mechanical energy of the polymer. The problem is usually considered in the 
following terms: at the first step, the thin rod which models the molecule is replaced by 
its centerline. With each point of the line (specified by its position along the contour C.) 
one associates a triad of unit vectors: the tangent to the line (t), the principal normal (n) 
which lies in the plane defined by the tangents at points ^ and ^ + dC,, and the binormal (b) 
which is orthogonal to both t and n. As one moves along the line, the triad rotates and this 
rotation is described by the Frenet-Serret equations in which the "rate" of rotation of each 
unit vector is determined by two parameters: the local curvature k, and the local torsion 
u (sometimes referred to as writhe)0. In order to relate this purely geometrical picture to 
the elastic response of real rods, one has to specify the physical properties of the rod in a 
stress-free (undeformed) reference state and to write the energy as a quadratic expansion in 
deviations from this state. In the classical theories of thin elastic rod^ one usually assumes 
that the reference state corresponds to a straight untwisted rod (with vanishing spontaneous 
curvature kq and spontaneous torsion uq) and the mechanical energy density is written as 
a sum of terms proportional to and u"^. The generalization to the case of non-vanishing 
spontaneous curvature and torsion is then done by requiring that the strain energy density 
per unit length U is minimized for k, = and u = uq which leads to the expression 



Here Ai and A2 are material parameters (products of elastic moduli and moments of inertia). 



qualitatively, the experimental data on torsionally constrained DNASJHfl. To account for the 
coupling between bending and twist observed in experiment, extra terms are conventionally 
added to the mechanical energy density, Eq. (1), by hand. 

The objective of this work is to derive an expression for the mechanical energy and obtain 
the equations which determine the mechanical equilibrium of a rod subjected to arbitrary 
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forces and moments. This is done using a new form of the displacement field, which accounts 
for both the deformation of the centerline and the rotation of the cross-section around this 
line (i.e., twist). Instead of using ad hoc assumptions about the form of the coupling between 
bending and twist, we will use standard methods of the theory of elasticity in order to derive 
the correct form of the coupling. 

In this work we will consider cylindrical rods with circular cross-sections. Although, 
at first sight, this case appears to be simpler than that of rods with asymmetric cross- 
sections, the reverse is true: while in the asymmetric case one can introduce a triad of 
vectors associated with the principal axes of inertia, which can rotate at a different rate than 
the Prenet triad, no such natural choice is possible in the symmetric case which therefore 
requires a more careful analysis. 

The exposition is organized as follows. Section 2 deals with geometry of deformation. 
The strain energy density of a rod is introduced in Section 3. Stress-strain relations are 
developed in Section 4. In Section 5, force and moment balance equations which describe 
the mechanical equilibrium of thin rods are derived. Several examples which illustrate the 
different aspects of the interaction between elongation, torsion and twist, are discussed in 
Section 6. Finally, in Section 7 we discuss the connection between our results and other 
theoretical and experimental works and outline directions for future research. 

2. Geometry of deformation 

A long chain is modeled as an elastic rod with length L and a circular cross-section S with 
radius a <^ L. Denote by ^ the arc-length of the centerline of the rod in the reference 
(stress-free) configuration. Let Ro(0 be the radius vector of the longitudinal axis and 
to(0 — dRo/d^ the unit tangent vector in the reference state. The unit normal vector no(^) 
and the unit binormal vector bo(C) are introduced by the conventional way. These vectors 
obey the Prenet-Serret equations with given spontaneous curvature ko(0 and torsion cuoiC)' 
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dto dno dho 

— = Kono, — = uoho - Koto, — = -a;ono. (2) 

Points of the rod refer to Lagrangian coordinates {^i}, where ^i, ^2 are Cartesian coordinates 
in the cross-sectional plane with unit vectors no and bo and ^3 = ^, 

rote, 6,0 = Ro(0 +6no(0 +6bo(0- (3) 

It follows from Eqs. (2) and (3) that the covariant base vectors in the reference configuration, 
gok = dro/d^k are given by 

goi = no, go2 = bo, go3 = (1 - ft:o6)to + ^^o(6bo - 6no)- (4) 

The position of the longitudinal axis of the rod in the actual (deformed) configuration is 
determined by the radius vector R(0- Following the conventional theories of rods, see, 
e.g.,il, we assume that the longitudinal axis is inextensible, which means that ^ remains 
the arc-length in the actual configuration (for attempts to account for the extensibility 
of the longitudinal axis, seelllt- The unit tangent vector in the actual configuration 

t = dH/dC, together with the unit normal vector n and the unit binormal vector b satisfy 
the Frenet-Serret equations 

dt dn , dh , . 

- = nn, ^=^h-nt, ^ = "^n. (5) 

For Kirchhoff ro dsi, the radius vector of an arbitrary point is represented as an expansion 
in the coordinates ^1 and 

r(ei, 6, = R(0 + ein(0 + e2b(0. (6) 

The functional form of Eq. (6) implies that any cross-section remains planar and perpendic- 
ular to the centerline of the rod, even in the actual deformed configuration. Furthermore, it 
also implies that any cross-section rotates rigidly with the longitudinal axis and therefore 
Eq. (6) does not allow for the possibility of a twist of the cross-section with respect to 
the centerline of the rod. Since the latter assumption has no physical basisi, we relax it by 
introducing a more general displacement field 



6,0= R(0 + (6 COS a - 6 sin Q;)n(0 + (6 sin a + ^ cos Q;)b(0, (7) 

where a{^) is the rotation angle around the centerhne of the rod. Prom here on we will refer 
to this rotation as "twist" and will reserve the terms "torsion" and "writhe" to describe the 
three-dimensional geometry of bending of the centerline of the rod. 
The covariant base vectors = dr/d^k are given by 

gi = cos an + sin ah, g2 — — sin an + cos ah, 
g3 = 1 — K,{C,i COS a — ^2 sin a) t 

'^{^'^~^) — (6 sina + ^2C0SQ;)n + (^1 cosa — ^2 sinQ;)b . (8) 

The contravariant base vectors g'^ are found from Eq. (8) and the equahty gj • g-' = Sf, 
where the dot stands for inner product and 5j is the Kronecker delta. Simple calculations 
result in 

g^ = cosan + sinab + Cit, g^ = — sinan + cosab + C2t, g^ = Cat, (9) 
where 

An = -{-77 + Co") (^1 sin a + {2 cos a), A = + a;)({i cos a - {2 sin a), 

At — 1 — k{^i cos a — ^2 sin a), Ci — — — {An cos a + A^sin a), 

At 

C2 = — (^nSino; - Abcosa), C3 = — . (10) 
At At 

One can proceed to calculate the energy of deformation using the displacement gradient 
either in the deformed, Vro, or in the reference, Vor, state. Both approaches result in 
the same expression for the mechanical energy. We will use the displacement gradient in 
the actual configuration Vro, because the corresponding strain tensor is connected with 
the stress tensor (always defined with respect to the coordinates in the deformed state) 
by conventional constitutive equations. It follows from Eqs. (4) and (9) that the tensor 
Vro = g'^gofe is given by 
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Vro = cos a(nno + bbo) + sin a(bno — nbo) 

+ {Ci - C3cuo6)tno + {C2 + Cso^oe^tbo + ^3(1 - K^ii)tio. (11) 

As a measure of deformation, the Almansi tensoril A = Vro ■ Vr,} is employed, where T 
stands for transpose. The tensor A is connected with the strain tensor e in the deformed 
state by the equahty e = |(I — A), where I is the unit tensor. It follows from Eq. (11) that 
the non-zero components of e are given by 

1/^ • ^ ./da \ 
ei3 = esi = --(4isma + 42Cosa)(— + u; - c^oj, 

1/^ ^ . .(da \ 

£23 = £32 = 2(^1 cos a - 42sma)(— + u; - t^oj, 

£33 = -«^(Ci cosa - ^2sina) + Ko^i, (12) 



where we kept only terms up to first order in ^1 and .^2- The neglect of second and higher 
order terms follows from the standard small local deformation assumption, which implies 
that all the length scales associated with bending, torsion and twist (e.g., radii of curvature) 
are much larger than the diameter of the rod. Note that this approximation is consistent 
with the form of the displacement field, Eqs. (3) and (7), where only terms up to linear 
order in the transverse coordinated ^1 and ^2 were kept. 

3. Strain energy density 

For a linear anisotropic elastic medium, the mechanical energy of elongation per unit volume 
in the deformed state is calculated as 

Mel = ^^14' (13) 

and the mechanical energy of shear is 

= ^2(e?3 + eL + 4 + 4), (14) 

where Ei and E2 are the appropriate elastic moduli. It follows from Eqs. (12) to (14) that 
the mechanical energy density 



M = Mel + Ush (15) 

is can be written as 

M = i|^i K\^icosa-^2smay + Kl^i^]+E2i^l + Q{^+uj-uoyy (16) 
The mechanical energy per unit length is given by 

which yields, upon integration 
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with 



'da 



Ai (^K^ - 2kko cos a + Kq) + ^2 (— + - cjq) (17) 



A, = EJ, A2 = 2E2I, 1=1 ildiidi^ = I ildiidi^. I iii2diidi2 = 0. 

Comparison of Eqs. (1) and (17) shows that the two expressions coincide in the absence of 
rotation of the cross-section with respect to the centerline (no twist, a = 0). In the general 
case, when a 7^ 0, Eq. (17) differs from Eq. (1) in several important ways: 

1. The torsion uj is replaced by u + da/d^. This correction has a simple intuitive meaning: 
the rotation of a point on the surface of a rod is the sum of the rotation in space of 
the centerline of the rod and of the twist of the cross-section about this centerline. 
Notice that this correction may always be present, independent of whether the rod 
has a non-vanishing spontaneous curvature (kq) and spontaneous torsion [ujq) or not. 
Such a correction was, in fact, proposed by previous investigator^. 

2. The term 2k,k,q is replaced by 2kko cos a, introducing a non-trivial coupling between the 
spontaneous and the actual curvatures of the rod, and the twist of its cross-section 
with respect to the centerline. Note that this term appears only when the rod has 
a non-vanishing spontaneous curvature and therefore while it has no effect on the 
elasticity of straight rods {kq = 0), it has a dramatic effect on the elasticity of helices 
and other curved (kq 7^ 0) rods. 



3. The usual expression for the energy, Eq. (1), is minimized when the curvature (k) and 
torsion (a;) recover their spontaneous values {kq and coq, respectively) in the stress-free 
reference state. Although this appears to be no longer true for our energy, Eq. (17), 
the difference stems from the fact that we have introduced a new independent variable 
(a) that describes the twist of the cross-section with respect to the centerline of the 
rod. In the absence of externally applied torques and tensile forces, minimizing the 
energy with respect to k, cu and a yields their values in the stress-free reference state, 
i.e., Ko, ujo and a — 0, respectively. 



4. Stress-strain relations 

Denote by a the Cauchy stress tensor and by cr*-' its contravariant components in the basis 
of the actual configuration. Substitution of Eqs. (13) to (15) into the equality 

du 



dsij 



results in 

(7^3 = 0-32 ^ ^^^^^ a -^2 cos a) + a; - a;o) , (18) 
(7^3 = Ei[-K{^iCOsa - ^2 sin a) + KqCi]- (19) 

Equation (19) does not take into account the inextensibility of the longitudinal axis. In 

order to enforce this constraint, we add an unknown parameter p (a Lagrange multiplier 
analogous to pressure for incompressible solids) to Eq. (19): 

0-33 = -pjf. El cos a - ^2 sin a) + KqCi] ■ (20) 

Since the unit normal to a cross-section of the rod coincides with t, the internal force (per 
unit area) f that acts on the cross-section of the rod is given by 

f = t • (7 = cr^^n + cr^^b + a^H. (21) 



It follows from Eq. (7) that the radius vector p from the center point of the cross-section 
(its intersection with the centerline) to an arbitrary point of the cross-section is 

p — (^1 cos a — ^2 sin Q;)n + (^i sin o; + ^2 cos a)h. (22) 

The moment (per unit area) p of the internal force with respect to the center point of the 
cross-section is defined as 

= p X f, (23) 

where x stands for vector product. Combining Eqs. (21) to (23) and using the equalities 

txn = b, nxb = t, bxt = n, (24) 

we find that 

p = [^(^1 sin a + ^2 cos a)n — {^i cos a — ^2 sin a)h a^^ 

+ |^(^i cos a — ^2 sin a)a'^^ — (^1 sin a + ^2 cos a)a^^ t. (25) 

The internal moment M is obtained by integrating p over the cross-section of the rod, 

M = pd^id^2 = M„n + M^h + M^t. (26) 

In principle, one could proceed in similar fashion and obtain the internal force 

F =F„n + F^b + Ftt (27) 

by integrating f over the cross-section of the rod. However, inspection of Eqs. (18)-(21) 
shows that since our expression for f is linear in the transverse coordinates ^1 and ^2, the 
integral over the cross-section vanishes. The source of the problem can be traced back to our 
choice of the displacement fields, Eqs. (3) and (7), where only linear terms in the transverse 
coordinates ^1 and ^2 were taken into account. Note, however, that even if we were to 
keep higher order terms in ^1 and ^2 in these equations, the unknown function F would be 
expressed in terms of new unknown functions (coefficients of quadratic contributions in ^1 
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and ^2 to the displacement fields). Instead, we will follow the standard approach and treat 
the vector F as an additional unknown that is found from the equilibrium equations (force 
and moment balance conditions). 

We now proceed to calculate the internal moment by substituting expressions (18) and 
(21) into Eqs. (25) and (26). Upon integration we obtain the constitutive relation between 
the parameters that characterize the deformation [k, uj and a) and the internal moment M 

M = AiKq sin an + Ai[k — kq cosa)b + A2(uj + — — uJojt. (28) 

Equation (28) is a new expression for the moment of internal forces which accounts for the 
twist of the cross-section with respect to the centerline of the rod. As expected, the internal 
moment vanishes in the stress-free reference state: k — kq, ou — cuq, a — 0. In the absence 
of twist, a = 0, Eq. (28) reduces to the conventional expression for Kirchhoff rods 

M = A(k - Ko)b + ^2(0; - a;o)t. (29) 



5. Equilibrium equations 

Consider an element of the rod bounded by two cross-sections with longitudinal coordinates 
^ and ^ -I- d^. Forces acting on this element consist of the internal force — F(^) apphed to 
the cross-section ^, the internal force F(^ -|- d^) apphed to the cross-section i + d^, and the 
external force qd^ proportional to the length of the element d^. Balancing the forces on the 
element yields 

F(e + dO - F(0 + q(Oc?e = 0. (30) 
Expanding the vector function 

F(e + dO = F„(e + dOnii + dO + F,{i + dOHi + dO + Ft(e + dOtii + dO 

into the Taylor series, using Eq. (5), and neglecting terms of second order in d^, we find 
that 
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+ KFt - ujFk)n + + ujFn)h+ - KFn)t 



Substitution of Eq. (31) into Eq. (30) results in the equilibrium equations 



dC 



(31) 



dR 



d^ 

m 

d^ 



dC 



+ ujFn + gft = 0, 



kF„ + = 0. 



(32) 
(33) 



where g^, and qt are the components of the external force per unit length, q = g^n + 
q^h + qtt. 

The moments acting on the element of the rod consist of the internal moment — M(^) 
applied to the cross-section ^, the internal moment M(^ + d^) applied to the cross-section 
^ -I- d^, the moments of internal forces — F(^) and F(^ -|- d^), and the external moment md^ 
proportional to the length d^, where m = m„n -|- mfcb -|- m^t is the external moment per 
unit length. To first order in d^, the moment of internal forces with respect to the center of 
the cross-section with coordinate ^ is 



[R(e + di) - R(o] X F(e + di) = m x F(e)de. 



The balance equation for the moments reads 



M(e + - M(e) + m X F(OrfC + m(C)rfe = 0. 



(34) 



It follows from Eqs. (24) and (26) that t x F = -F^n + F„b. By analogy with Eq. (31), 
one can write 

^dMt 



M(e + dO - M(0 



[-7^ + i^Mt - ouMb) n + (— ^ + uMn) h 



Substitution of these expressions into Eq. (34) results in the equations 



dC 



dMn 

dC 
dMt 

d^ 



+ KMt - cuMb - Fft + mn = 0, 
- kM^ + mt = 0. 



dMb 
~d^ 



+ ujMn + F„ + mfe = 0, 



(35) 
(36) 



Given the vectors M and m, Eqs. (35) can be used to determine the forces F^ and F^. 
Eliminating the unknown functions F„ and F^ from Eqs. (32) and (35), we obtain 
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^ + i^(^^ + ujMn + mb^+qt = 0, (37) 

^ (^-^ + i^Mt - uMh + rrinj - ^[-^ + + m^) + = 0. (38) 

Equations (36) to (38) together with constitutive relation (28) are a set of four nonhnear dif- 
ferential equations which determine the four unknown functions Ft, a, k and uj. The neglect 
of a (that is the use of conventional formula (6) instead of Eq. (7) for the displacement field 
r) is acceptable only when special restrictions are imposed on external forces and moments. 
In the general case, this simplification is not correct, and Eq. (7) should be employed for 
the analysis of deformations. 

6. Examples 

A. Twist of a closed loop 

Consider a rod whose stress-free shape is a planar circular loop with radius oq, under the 
action of a constant twisting moment m^. It is assumed that the moments m„ and m^, as 
well as the forces and qt vanish. The solution of Eqs. (36) to (38) reads 

K — Kq — Oq ^, (39) 

2 

o; = cuo = 0, Ft = 0, a = arcsin (40) 

According to these equalities, any cross-section of the rod twists around its centerline by a 
constant angle a. This solution is not described by the Kirchhoff theory of thin rods. It 
exists as long as the moment rrit satisfies the condition |mt| < Aio^'^. If the latter restriction 
is not fulfilled, the planar shape of the loop becomes unstable. 
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B. Torsion of a disconnected ring 



We analyze the deformation of a disconnected ring (no contact between the points ^ = 
and C, = L). The end ^ = is fixed, and a torque T is apphed to the free end ^ = L. 
The centerhne of the rod in the stress-free reference state describes a planar circle with 
radius = and no spontaneous torsion, uq = 0. Similar problems were recently studied 
m and their solutions were applied to the analysis of kink transitions in short DNA 

rings. 

We assume that in the deformed state the centerline becomes a non-planar curve whose 
radius of curvature remains unchanged, see Eq. (39). For simplicity, we confine ourselves to 
small displacements, and neglect terms of order in the constitutive equations (28). This 
yields 

M„ = AiKoa, Mfe = 0, Mt = + uj). (41) 

Substitution of these expressions into the equilibrium equations (36) to (38) implies that 
the longitudinal force Ft vanishes, whereas the functions a and uj obey the equations 



'd'^a duj\ . o d 



da . /da 



. /u a auj\ .9 u , uu; , /ua \ , 

^(^ + -)-A.^a = 0, ^^A,- + A{- + u) =0. (42) 



d^ ^^d^ 



It follows from the second equality in Eq. (42) that 

{Ai + A2)^ + A2UJ = c, (43) 

where c is a constant to be found. Excluding u from Eqs. (42) and (43), we obtain 
d'^a 2 



+ nla = 0. (44) 



The solution of Eq. (44) is given by 

a = Ci sin + C2 cos kq^, (45) 

where ci and C2 are arbitrary constants. Substitution of Eqs. (43) and (45) into the boundary 
conditions at the clamped end C, = 
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a(0) = 0, uj{0) = 
implies that 

a = CisinKo^, = — — -kqCi{1 — coskq^). (46) 

Equating the moment at the end ^ = L to the external torque T and using Eqs. (41) 
and (46), we obtain 

T 

which results in the formulas 
T 



a 



A2K0 ^ " ~^~34|~~^(^ ~ cos(ftoO)- (47) 

Equations (41) and (47) provide an explicit solution to the torque problem, which cannot 
be obtained in the framework of the Kirchhoff theory of rods. When the radius of the ring 
tends to infinity, i.e. for a prismatic rod, Eq. (47) implies that 

« = ^e, ^ = 0. (48) 

In this limit, the solution (48) coincides with the classical displacement field for the twist of 
a circular cylinder^. 



C. Helix under tension and torque 

A helix-shaped rod whose stress-free reference state is characterized by spontaneous curva- 
ture Ho and torsion uo, is deformed by tensile forces P and torques T applied to its ends. All 
other forces q and moments m are assumed to vanish. We introduce Cartesian coordinates 
{xk} with unit vectors and describe the configuration of the centerline of the rod in the 
stress-free reference state by the vector 

Ro = ao cos ei + oq sin 62 H — . 63. (49) 

/al + bl ^al + hi ^al + hi 

14 



The parameters oq and 60 are expressed in terms of the spontaneous curvature kq and torsion 
cuq by the formulas 



1. Fixed force and torque on ends 

Consider a rod whose centerhne describes one complete turn of a helix (the angle between 
tangent vectors at the two ends of the undeformed rod equals 27r). The contour length of 
the rod is 

I ^ 27r(i4 + u^yK (51) 
We assume the following boundary conditions at the ends of the rod: 
M„(0) = M„(0 = 0, Mh{0) = Mb{l) = 0, 

M,{0) = M,{l)=T, Fi(0) = Fi(/) = P. (52) 

Equations (52) imply that the torque T and the tensile force P are the only external loads 
applied to the segment. Assuming the parameters P and T to be rather small and neglecting 
the deviation of torsion from its value in the stress-free state, we look for a solution of the 
equilibrium equations in the form 

a — Aa, K — Kq + Ak, cu — cuq, (53) 

where Aa is small compared to unity, and Ak is small compared to kq- 

Neglecting terms of the second order in the perturbations of twist angle and curvature 
{Aa and Ak, respectively), we find from Eq. (28) that 

Mn = A^KoAa, Mb = A^Ak, Mt = ^2^. (54) 

We substitute expressions (53) and (54) into Eqs. (36) to (38), neglect terms of the second 
order in Aa and Ak, and arrive at the equations 
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dMt 

— KoMn = 0, 

dFt ^dMb , X 

d'^Mb ^ dMn T. r. 

-—f + 2ujo—^ + KQUJoMt - ujlMb - KoFt = 0, 

d'^Mn dMt dMb 
d^^ di d^ 



(55) 
(56) 
(57) 
(58) 



where the longitudinal force is assumed to be small as well. It follows from Eqs. (55) and 
(58) that 



dMu 



d^ 2ujq 



+ {kI - ul)Mr, 



(59) 



Substitution of Eq. (59) into Eq. (56) results in 



dFt , 



+ 



d^ 2ujq 



- + {4+ ul)M^ 



= 0. 



(60) 



Equations (55), (57) and (59) imply that 

dFt _ 1 (d'Mb , ^ d'^Mn , dMt ^dMb 

d^Mn / 2 I o 2\d'^^n , 2/ 2 , 2\i\/r 
+ {Kq + 2ujQ)—— + UJo{KQ + UjQ)Mn 



2KoUJo 



(61) 



Excluding the function Ft from Eqs. (60) and (61), we obtain a closed equation for the 
internal moment M^ 



The solution of Eq. (62) reads 

M„ = (ci + c'lO sinynl + ujU) + (c2 + 40 cos(^«g + a;^^), 



(62) 



(63) 



where c^, are constants to be found. It follows from the boundary conditions (52) for the 
function M„ and Eq. (63) that 



C2 = 4 = 0. 



(64) 



16 



Integrating Eq. (59) from to I and using boundary conditions (52) for the function M^, 
we obtain 



Substitution of expressions (63) and (64) into this equahty results in 

c[ = 0. (65) 
Combining Eqs. (54) and (63) to (65), we find that 

A«(0 = -^smUKl + coU), (66) 

Note that although the twist angle vanishes at the ends and in the middle of the rod 
(Aa(0) = Aa(/) = Aa(//2) = 0), it does not vanish elsewhere. Differentiating Eq. (66) and 
using Eq. (54) and the boundary conditions (52) for M^, we arrive at the equality 

ci = , (67) 



Substitution of Eqs. (64), (65) and (67) into Eqs. (54), (63) and (66) implies that 



Mn^ ^^L=sin(J'^^li), M,^Tcos(yjKl-ruli). (68) 



It follows from Eqs. (59) and (68) that 



Integrating this equality with the boundary conditions (52) and substituting in Eq. (54) 
yields 



A«(0 = 



1 



(69) 



Note that the sign of Ak vanishes at the ends of the rod; inside it, its sign is opposite to 
that of the torque T (positive torque means overtwisting) . Substitution of Eqs. (68) and 
(69) into Eq. (57) gives the internal tensile force 
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F, = Tc.o[(l + ^(l + 



'0^ 



l))cOs(yK2T^^) 



(70) 



It follows from Eq. (70) that our solution ^ and q;(^) under boundary conditions (52) 
is valid if the tensile force P and the torque T applied to the ends of the rod satisfy the 
relation 



Equations (68) to (70) provide an explicit solution to the problem of combined tension and 
torsion of a helical segment. The main results are as follows: 

1. The application of positive torque T at the ends (overtwist) leads to axial compression 
of the helix which is maximal at the center and vanishes at the ends of the rod; 

2. The ratio of the tensile force P and the torque T is independent of the initial curvature 

(and, therefore, of the length of the rod) and depends only on the initial torsion uoq 
and the ratio of elastic moduli A1/A2 = Ei/{2E2); 

3. The force P is proportional to the torque T. This result is markedly different from 
that obtained for a similar deformation of a circular incompressible cylinder, where P 
can be shown to be proportional to (the Poynting effectS). 




Note that our solution corresponds to a helical rod (with constant and ujq) which, 
upon application of external forces and torques, is deformed into a new, non-helical shape. 
It is natural to ask under which boundary conditions a helix will deform into another helix 
(with constant k, and u), and derive the corresponding force-elongation relation. This is 
done in the following. 

2. Elongation and winding of a helix 

Consider a helix made of an arbitrary number of repetitive units Co {Cq is the smallest 
segment for which the angle between tangent vectors at its ends is 2tt) such that its stress- 
free reference state is characterized by the parameters kq and ujq. We allow deformations 




(71) 
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that satisfy the following conditions: (i) the rod becomes a helix with constant curvature k 
and constant torsion cu, and (ii) the twist a vanishes. Under the action of combined tensile 
force P and torque T, any repetitive unit Cq of the helix in the reference state is transformed 
into an element with the angle between tangent vectors at the ends 27r(l + </?), where the 
angle 2nip can take positive or negative values. The radius vector of the centerhne of the 
rod in the deformed state can be written in the form 

R = a cos{Si)ei + a sm{Si)e2 + Si^e^, (72) 

where a, S and 5*1 arc constants which will be calculated in the following. Differentiating 
Eq. (72) with respect to ^ and bearing in mind that |t| = 1, we obtain 

a^S^ + = 1. (73) 

According to the definition of ip, 

SI = 27r(l + (p). (74) 

The projected distances (along the xs-axis) between the ends of the repetitive unit in the 
reference and deformed states are 

Ho = 2nbo, n = Sil, (75) 

respectively. The axial elongation r] is defined as the ratio of these distances, 
n Sil I 



Simple calculations result in the formulas 

= a,5^ uj = Sy/l - a?S'^. (77) 

It follows from Eq. (72) that the projection of the force Ft on the axis x^, is Fjt • 63 = FfSi. 
Equating this expression to the tensile force P and using Eq. (73), we arrive at the relation 

Ft = , ^ (78) 
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which means that Ft is independent of ^. Equation (28) imphes that components of the 
moment M are independent of ^ as well, 

Mn = 0, Mfe = Ai{k - fto), Mt = A2iuj - uo). (79) 
The only equihbrium equation reads 
uj{KMt - tuMb) = KFt. 

Substitution of expressions (78) and (79) into this equality yields 



LU 



A2k{uj - ujq) - Aiuj{k - Ko) = (80) 



Excluding the parameters a, S, k and u from Eqs. (73), (74), (76), (77) and (80), we express 
the tensile force P in terms of the axial elongation of the helix rj: 



-To 



where 



;i + VP)r^ - 1 - Ar^{{l + ^) - , )] , (81) 



Comparing Eq. (52) with Eqs. (78) and (79) we find that the only difference between the 
boundary conditions for the two problems is that in the former case we have neglected the 
moment Mf,. The fact that a minor change of boundary conditions can drastically change 
the character of deformation is quite remarkable and indicates that these conditions should 
be chosen with care. 

Since various variants of the theory of elastic rods were applied to interpret the experi- 
mental force-elongation curves for stretched DNA molecules at large deformations, we will 
present plots of some of the results of this section and comment on their qualitative features. 
The graph Pq = Po{t]) for extension without torsion, ^9 = 0, is plotted in Figure 1. In the 
calculation we used A = 0.67, in agreement with conventional data on DNA Ai = 50 nm, 
A2 = 75 nmEl, where A^ = Ak/ (kBT), is Boltzmann's constant and T is temperature. 
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No detailed comparison with experiment is attempted here, but Figure 1 captures rather 
well the qualitative features of the experimental data for DNA moleculesHH. 

In order to check whether our theory captures the qualitative features of experimental 
data on the elasticity of supercoiled DNA, the dependence i] = rj{ip) is depicted in Figure 2 for 
various tensile forces Pq. This figure also shows qualitative agreement with observations on 
the DNA chainJnii^: for small tensile forces, there is pronounced asymmetry with regard to 
the sign of (p, but the r] = ri{(f) curve becomes nearly flat at large tensile forces. Throughout 
the parameter range, the elongation decreases nearly linearly with degree of supercoiling. 
All these features were observed experimentally and were interpreted as a proof for the 
existence of a new type of twist-stretch coupling0'Q. Note, however, that in the analysis 
that led to Figure 2 we assumed that the deformation of the helix takes place with no twist 
of its cross-section around the centerline of the rod {a = 0). Therefore, our solution can 
be derived using the standard theory of elastic rods, based on the elastic energy of Eq. 
(1), in which no such coupling appears. Inspection of the derivation of Eq. (81) leads to 
the conclusion that the strong dependence of elongation on the degree of supercoiling has 
a simple physical meaning: when an inextensible helical rod is subjected to torque that 
produces supercoiling, each new turn has non-vanishing projection on the Xi — X2 plane and 
the projection of the deformed helix on the X3 axis (i.e., its elongation) decreases as the 
result. This can be fully described by Eq. (1) and does not require the introduction of new 
coupling into the mechanical energy of elastic rods. 

7. Concluding remarks 

In this work we have extended the theory of elasticity of thin inextensible rods beyond that 
of three-dimensional space curves which can be completely described by local curvature 
K and geometric torsion u. We have shown that in order to describe the displacement of 
a point in a rod of arbitrarily small but non-vanishing thickness, one has to account for 
deformations that produce a rotation of the cross-section of the rod about its centerline. 
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The modified displacement field was then used to calculate the strain tensor. The resulting 
expression for the mechanical energy of rods with non-vanishing spontaneous curvature 
contains a new coupling term between the curvature of the rod and the twist of its cross- 
section with respect to the centerline, which does not appear in any of the previous theories. 
We derived the complete set of non-linear differential equations which describe the conditions 
of mechanical equilibrium and which can be solved for the parameters of deformation n, uj 
and a for arbitrary external forces and moments acting on the rod. In order to illustrate the 
physical consequences of our theory, we proceeded to analyze several illustrative examples. 
In particular, we have analyzed the deformation of a helical rod subjected to a combination of 
tension and torque and showed that the theory captures the qualitative features of the recent 
observations on the connection between supercoiling and elongation of strongly stretched 
DNA molecules. 

Note that we have described the deformation of thin rods by three independent functions 
a, . and ... This is reminiscent of the conventional approach! where the deformation is 
described in terms of the three components of the, so called, "twist" vector, ki, K2 and K3. 
Although this was not mentioned by the above authors, such an approach goes beyond the 
purely geometric description of an elastic line in which only two functions are necessary!^ 
and describes a line with some "internal structure". With each point of this line one can 
associate a "physical" triad of vectors that differs, in general, from the "geometric" (Frenet) 
triad. While the two triads have one common vector (the tangent to the line), the other two 
pairs of vectors rotate at different rates as one moves along the line contour and therefore the 
rotation of the physical triad can not be completely described by the two Frenet parameters 
K and UJ. It is important to realize that the introduction of a physical triad is necessary 
whenever some asymmetry of the cross-section, either geometric or physical^, is present. 
However, even though the procedure is not unique, one may also introduce the physical triad 
by hand even for a rod with a circular cross-section. For example, we may draw a line on 
the surface of a rod which describes the intersection of the normal vector with this surface. 
When the rod is deformed, the deformation of this line will, in general, be different from 
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that of the centerhne. We can now connect the corresponding points of the two hnes (having 
the same contour parameter ^) and define the resulting vector as one of the vectors of the 
physical triad. The remaining vector is then defined as the normal to the plane formed by 
the above vector and the tangent to the centerline. This procedure is completely equivalent 
to what we have done here, by introducing the rotation a{^) and explains the appearance of 
an a— dependent term (k kq cos a) in the expression for the mechanical energy, that couples 
the curvatures in the stress-free and the deformed states of the rod. Note that while for 
rods with asymmetric cross-sections, three independent parameters are needed in order to 
characterize the stress-free reference state, only two such parameters (e.g., kq and uq) are 
necessary in the degenerate case of rods with circular cross-sections. 

There are several possible directions in which the work presented here can be extended. 
For example, throughout this work we assumed that the conditions of mechanical equilibrium 
can be satisfied and considered only stable configurations of the deformed rods. However, the 
introduction of a new type of deformations is expected to have a profound effect on various 
instabilities (e.g., buckling under torsion and twist, plectoneme formation, etc.) and we are 
now studying these questions. Another direction for future research involves the extension 
of the present, purely mechanical, analysis to include the effects of thermal fluctuations. 
This leads naturally to a new class of physical models for rigid biopolymers and protein 
assemblies which can account for the spontaneous curvature^ of these objects. 
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FIGURES 

Fig. 1. The dimensionless tensile force Pq versus the axial elongation rj for a helix with A = 0.67 
and if = 0. Curve 1: A = 0.5; curve 2: A = 0.6; curve 3: A = 0.7 

Fig. 2. The axial elongation i] versus the overtwist (p for a helix with A = 0.67 and A = 0.6. 
Curve 1: Pq = 0.1; curve 2: Pq = 1.0; curve 3: Pq = 10.0 



26 




Fig. 1. The dimensionless tensile force Pq versus the axial elongation t] for a helix with A = 0.67 
and = 0. Curve 1: A = 0.5; curve 2: A = 0.6; curve 3: A = 0.7 
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Fig. 2. The axial elongation r] versus the overtwist (p for a helix with A = 0.67 and A = 0.6. 
Curve 1: Pq = 0.1; curve 2: Pq = 1.0; curve 3: Pq = 10.0 
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